Magnetic Phases of Electron-Doped Manganites 
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We study the anisotropic magnetic structures exhibited by electron-doped manganites using a 
model which incorporates the double-exchange between orbitally degenerate Cg electrons and the 
super-exchange between t2g electrons with realistic values of the Hund's couphng(Jj^), the super- 
exchange coupling( Jaf), and the bandwidth(W). We look at the relative stabilities of the G, C and 
A type antiferromagnetic phases. In particular we find that the G-phase is stable for low electron 
doping as seen in experiments. We find good agreement with the experimentally observed magnetic 
phase diagrams of electron-doped manganites [x > 0.5) such as Ndi^ajSr^^MnOs, Pri-^^Sra^MiiOa, 
and Smi_a;Caa;Mn03. We can also explain the experimentally observed orbital structures of the C 
and A phases. We also extend our calculation for electron-doped bilayer manganites of the form 
R2-2a: Ai+2i:Mn207 and predict that the C-phase will be absent in these systems due to their reduced 
dimensionality. 

PACS : 75.30.Vn, 75.30.Gw, 75.30.Et, 71.27.-t-a 



I. INTRODUCTION 

Recently there has been a great upsurge of interest in 
doped manganites exhibiting colossal magnetoresistance. 
Most of the studies on manganites ( Ri-xAxMnO-i ; R 
= La, Nd, Pr or Sm and A — Sr, Ca, Ba or Pb) focus 
on the transition from a ferromagnetic metal to a para- 
magnetic insulator in the doping regime 0.15 < a; < 0.4 
[Q-p| which can be considered as doping the quarter filled 
eg band of RMnOa with holes. More recently, electron- 
doped manganites [Q, namely systems wih 0.5 < x < 1, 
have begun to be explored. These seem to be different 
and quite interesting in their own way with a variety of 
anisotropic magnetic phases and with no evidence of par- 
ticle - hole symmetry. Such systems are experimentally 
seen to exhibit A-type, C-type and G-type antiferro- 
magnetism spectacular transitions between these 

phases under an applied magnetic field |^ as well as pos- 
sible phase separation. In contrast with the hole-doped 
systems, there have been very few attempts to under- 
stand the electron-doped systems theoretically. We show 
here that the rich magnetic phase diagram as well as 
their orbital structure can be understood in terms of a 
microscopic model which takes into account the large but 
finite Hund's rule double- exchange(DF,) coupling, effects 
of orbital degeneracy and the super- exchange (SE) cou- 
pling between t2g spins within a band picture. 

We study this model in detail for realistic values of the 
Hund's coupling, the super-exchange coupling and the 
bandwidth. We present the phase diagram as a function 
of doping in reduced units of Jn/t and JAp/t. From 
the phase diagram we deduce that the key interaction 
responsible for the stability of the G-phase near a; = 1.0 
is the super-exchange interaction. We also find that the 
A-phase near x 0.5 is very sensitive to the variation 
of the super-exchange interaction. We obtain a G-type 



phase for 0.85 < x< 1, a C-type phase for 0.6 < x< 0.85 
and an A-type phase for 0.5<a;<0.6 for values of Jh 
(Hund's coupling), Jaf ( super-exchange coupling ) and 
t (hopping parameter) that are in agreement with den- 
sity functional calculations Thus we find that a 
finite value of Jh leads to a magnetic phase diagram in 
good agreement with experiments for a number of doped 
manganites Ri-xAxMnO^ for x > 0.5. The model also 
throws light on the nature of orbital occupation of the 
electronic degrees of freedom which will lead to the exper- 
imentally observed orbital ordering We extend the 
mean field theory to incorporate the CE-phase at a; = 0.5 
and find that it is stabilized over a wide range of values 
of JHSo/t and JafSq/I. We also use our model to make 
a number of predictions regarding the magnetic phases of 
electron-doped bilayer systems R2-22;Ai+2a;Mn207 [ p^ . 
Specifically, we point out that the C-phase will be absent 
in the electron doped bilayer manganites due to reduced 
dimensionality. 

We start by briefly describing the experimental situa- 
tion. In Section II we present our model. In Sections HI 
we present our results on the magnetic phase diagram 
of the manganites. Section IV deals with the nature of 
the G-phase for low electron doping and of canting of 
spins. In Section V we explain the orbital structures ob- 
served in the C and A phases. Section VI deals with the 
phase diagram of the electron-doped bilayer manganites. 
Finally we make a comparison of our results with the 
earlier works and point out the shortcomings of various 
approaches including ours. 

The conventional single-band double exchange model 
predicts a phase diagram symmetric about x = 0.5. 
However, the behaviour of the observed ground state 
magnetic properties does not agree this simple picture. 
Experiments show a remarkable asymmetry with re- 
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gard to the magnetic properties of the system. In 
Lai-xSr^MnOs an A - type antiferromagnetic ground 
state is seen for 0.52 < a; < 0.58, above which it be- 
comes a C-type antiferromagnet. In Ndi-^Sr^MnOs an 
A- type antiferromagnetic phase extends from x ^ 0.5 
to X ^ 0.62 and a C-type phase is seen till x ^ 0.8 
In Pri-xSr^MnOs the A-type antiferromagnctism 
is seen from x ~ 0.48 upto x ^ 0.6 and the C- 
antiferromagnetism upto x ^ 0.9 The end compound 
AMnOs is a G-type antiferromagnet and this state ex- 
tends, in general, upto x ~ 0.90 |T^. In particular, re- 
cent experiments on Smi_2;Caa;Mn03 suggests that the 
G-phase, albeit with ferromagnetic clusters embedded in 
them, survives upto a doping concentration of a; = 0.88 
1^ . Though a picture based on band-structure will not 
be appropriate in such a case, we believe that the nature 
of the background magnetic phase can still be captured 
since the dominant energy here is the antiferromagnetic 
energy resulting from the super-exchange. 



II. MODEL 

If one starts from AMnOs and increase the doping, 
the doped electrons go into empty Cg levels doubly de- 
generate in the absence of Jahn- Teller splitting. As noted 
in the double-exchange between these degener- 

ate Cg levels along with the super-exchange between t2g 
core spins lead to a qualitatively different phase diagram 
which is highly asymmetric about x = 0.5. However, the 
resulting T = phase diagram they obtained while asym- 
metric is in disagreement with experiments on several 
counts (see below). This motivates a detailed study of 
an orbitally degenerate double-exchange (DE) and super- 
exchange(SE) model for 0.5 < a; < 1 with realistic values 
of parameters Jh,Jaf and bandwidth W [^sj. The ef- 
fective Hamiltonian describing the low energy properties 
of the system is 
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and (1^2 _y2 orbitals 



Here a and P denote the d^z^-r 
respectively. Si is the t2g spin at site i, Jh is the Hund's 
coupling and Jaf the super-exchange between t2g spins 
at nearest - neighbour sites i and j, and fi stands for the 
spin degree of freedom of the itinerant electrons. The 
hopping matrix elements are determined by the symme- 
try of eg orbitals . 

We treat the spin subsystem quasiclassically. As- 
suming a homogeneous ground state we take Si — 
So cos(Q.ri). where Q = (0,0,0) for the ferromagnetic 
phase, Q — (tt, tt, tt) for the G-type antiferromagnetic 



phase, Q = (tt, tt, 0) for the C-type antiferromagnetic 
phase, and Q — (0,0, tt) for the A-type antiferromag- 
netic phase. Canting can be included by assuming 
Sj — So{sm 9i, sin 9i, cos 9i) with 9i taking values be- 
tween and TT. This is discussed late below [ pT[ . 

Under these assumptions the electronic part of the 
Hamiltonian reduces to 
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-t{coskx + cos ky) — -tcoskz, 
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The super-exchange contribution to the Hamiltonian is 
given by 



JafSq 



{2 cos 9xy + cos9z) 



(4) 



with 0j.y — 9z — for ferromagnetic, 9xy = 9^ = tt for 
the G -type antiferromagnetic, 9xy = tt and 9z — for 
the C -type antiferromagnetic, and 9xy = and 9z — 
TT for the A -type antiferromagnetic phases. Here 9xy 
is the angle between nearest neighbour spins in the x- 
y plane, and 9z is the angle between nearest neighbour 
spins in the z direction. Inclusion of canting by assum- 
ing Si = 5*0 (sin sin cos 0i) will connect different spin 
species at the same site. These contibutions come from 
the (Tx and ay terms in the DE part of the Hamiltonian 
which are absent when canting is absent. Thus the DE 
part of the Hamiltonian becomes 



Hde = -JhSo E ^3 (4a 



JHSoy^,sm9j{c]^^Cjai + c]„^Cj„-r) 



JhSq i sm9j{c]^^Cjai - c]^^Cj„t) (5) 



We have neglected the correlation term in the Hamil- 
tonian UJ2ia'^haT^i.al + ^ Si.o-,!/ "■^i,i^"»,2,<t and the 

Jahn - Teller contribution gJ2i,a,i3.a4.a,aQi'^(^i,'y,P '^'^^^ 
Q describing the local distortion which lifts the degener- 
acy p8| . We neglect the correlation term because of the 
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low electron doping regime we are interested in. (a; = 0.5 
refers to a filling of 0.125 in our model and the filling 
ranges from to 0.125). For the same reason the inter- 
site Coulomb correlations, which may be necessary for 
the stability of charge ordered phase around x — 0.5, 
are also neglected. A cooperative Jahn- Teller effect can 
drastically change the magnetic ground state ||l^]. How- 
ever, since the carrier concentration is very small, so is 
the effective number of Jahn - Teller centres and hence 
we do not expect any qualitative change in the magnetic 
phase diagram though both may be required alongwith 
the breathing mode distortions induced by holes to ex- 
plain the CE - type charge ordered phase at x — 0.5 
Doping-induced disorder can have two effects. Firstly 
substitutional disorder may localize Cg electrons. How- 
ever as long as the localization length is more than the 
inter atomic spacing, the hopping to nearest neighbour 
sites will split the energy levels into bonding and anti- 
bonding orbitals with electrons occupying the bonding 
orbitals. This process is naturally taken care of in our 
model. Secondly the presence of a magnetic rare-earth 
ion can have coupling with the magnetic Mrv''^ ion and 
thus leading to change in the Mn-RE coupling as doping 
varies. However in most of the manganites, the RE ion in 
general is non-magnetic {eg. La) except, say in Pr. How- 
ever studies on Pr-Sr system around x = 0.37 pl[ have 
shown that Mn-Pr coupling plays no role in the magnetic 
propertis. Hence we also do not expect substitutional dis- 
order to play a role in determining the magnetic phases 
though, as argued in |13 1^, it might play a role in the 
transport properties of these compounds. We obtain the 
magnetic phase diagram by minimizing the total energy 
i/e; + HsA as a function of filling by fixing the chemical 
potential. 

We present the magnetic phase diagram for both man- 
ganites and bilayer manganites as a function of Jh /t and 
JAp/t. From density functional studies we estimate 
t = 0.15eV, JhS = 0.75eV and JafS^ = 8meV. Hence, 
we choose JnSo/t = 5 and JAFS^/t — 0.053. This value 
of JAFS^/t also leads to the correct mean field T/v for 
the end compound CaMnOa. Hence we use the phase 
diagram corresponding to these values for making com- 
parison with experiments |22(| . 



III. PHASE DIAGRAM OF THE ELECTRON 
DOPED MANGANITES 



The X — \ limit corresponds to empty Cg orbitals. The 
only contribution to the Hamiltonian comes from the SE 
interaction which is isotropic and hence leads to the G- 
phase at a; = 1. At low electron doping, however, the 
SE still wins over the Hund's coupling and leads to the 
G-phase. 



Doped electrons go into states with minimal energy 
corresponding to the F point at k = 0. (This is a conse- 
quence of the band-picture we, as well as other workers, 
use. The possibility of the doped electron forming ferro- 
magnetic clusters is mentioned later, in which case the 
band-picture will break down.) We first assume that all 
doped charges go into the state with k = and neglect for 
the moment the effects due to finite filling of the bands. 
(Strictly speaking this is only the case for very small dop- 
ing). Assuming uncanted states the energies for various 
magnetic states at k = are 



Eg = ~iJAFSll2t - y^l6+{JHSo/2ty, (6) 

Ea - JAFS^/2t -Ay- JHSoy/2t, (7) 

Ec = -JAFSl/2t - 8y/3 - y^l6/9 + {JHSo/2tf, (8) 

Ef = 3JAFSl/2t -4y- JHSoy/2t. (9) 

Here y is the actual electron filling in the two band model 
and is related to -4y + 1. For JnSo/t = 5 and 

JAFSfJt — 0.053 we find that the G-phase is stable upto 
X = 0.76 beyond which the A-phase becomes stable. In 
Fig(l) we present the phase diagram assuming the elec- 
trons go into the F point and there is no canting of core 
spins. We plot the phase diagram as a function of doping 
and JAFS'^/t for a fixed value of JHSo/t = 5 and doping 
and JHSo/t for a fixed value of JafSq/I = 0.053. 

The effects due to finite band-filling will alter these 
values and the numerically obtained values can be read 
from Fig(2) and Fig(3). This also leads to the physi- 
cally expected result that the doping region over which 
the G-phase stabilizes grows with JAp/t- As electron 
doping increases the kinetic energy starts dominating 
over the SE contribution leading to increased spin align- 
ment. This happens because kinetic energy is an in- 
creasing function of doping and for small doping it is 
proportional to the electron filling. However, a three- 
dimensional antiferromagnctic spin alignment does not 
allow for the motion of electrons. So to take advantage 
of the kinetic energy gain the mobile electrons polarize 
the spins along chains, planes and finally in all three 
directions successively. The reduction in the DE energy 
due to such alignment is overcome by the gain in kinetic 
energy beyond some doping value (for given values of 
Jn/t and JAF/t) and this point defines the G-C phase 
boundary. Moreover, as we will see in the next section, 
the C-phase has orbital ordering of c?z2-type and the 
A-phase has orbital ordering of dj.2_y2-type. Thus the 
interplay of the spin alignment along chains or planes 
and the corresponding orbital order also leads to change 



of the ^effective hopping parameters', t^ and t 



xy, 



the 



and X — y directions. In general, this leads to the system 
transforming from one-dimensional, to two-dimensional 
and finally three dimensional ferromagnetic structures 
with increasing doping. Thus the competition between 
effective kinetic energy (determined by Jh and band 
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filling) and super-exchange leads to transitions G-C-A- 
F(with number of antiferromagnetic bonds 6, 4, 2 and 
respectively) as the doping is varied for a given J///t. 

In Fig(2) we present the results for JhSq/I = 5. For 
JafS^/I = 0.053 we find that the system has a stable 
ferromagnetic ground state upto x = 0.47, the A-phase 
is favoured for x < 0.57, the C-phase upto 0.85. The 
G-phase becomes the stable phase for 0.85 < a; < 1. We 
also find that the A - phase near x = 0.5 is stable only for 
a limited range of JApSQ/t. The overall phase diagram is 
in excellent agreement with the experimentally observed 
phase diagram of NdSr, PrSr and SmCa systems. 

In Fig(3) we present the results for JafS'^/I = 0.053. 
We find that the A-phase, stable near x = 0.5 for smaller 
values of JuSo/t, gets pushed to the right making the 
ferromagnetic state stable near x = 0.5 for large values 
of Jh /t- However, in contrast to the earlier case, the A - 
phase is stable over a wide range of values of JhSq ft. We 
conclude that the A - phase near x = 0.5 is very sensitive 
to the variation of JafSq/I and rather less sensitive to 
the variation of JHSo/t. 

At X = 0.5 most of the manganites have a 
charge/orbital ordered ground state with the magnetic 
phase being the CE-type antiferromagnet. Our mean 
field theory can be extended to include the uncantcd CE 
phase by assuming ~ So/2 (cos((5i.ri) — cos{Q2.ri) + 
cos(Q3.ri) + cos((34-n)) with Qi = (7r,0, tt), Q2 = 
(0,7r,7r), Q3 = (7r/2, 37r/2, tt) and Q4 = (37r/2, 7r/2, tt). 
We present the phase diagram at x = 0.5 as a func- 
tion of JuSolt and JafSq/I in Fig(4). We find that 
at JHSo/t = 5 and JafSIH = 0.053, the CE phase 
stabilizes over other phases. In fact, the GE phase is 
stabilized over a wide region of the phase diagram at 
X = 0.5. This may explain why most of the manganites 
at X = 0.5 have the CE phase as the magnetic ground 
state. However, it is to be noted that the CE phase 
we obtain is not charge/orbital ordered. Other interac- 
tion such as the strong Coulomb repulsion or coupling 
of the lattice degrees of freedom to the eg electrons may 
be needed to make this phase charge/orbital ordered. 
There are contrasting views regarding the origin of the 
charge/orbital ordered CE phase and the precise role of 
JT and Coulomb effects is still not clear. Strong on-site 
Coulomb correlations within a two-band model seem to 
stabilize the CO phase at x — 0.5 |23|. It can also be 
thought of as emerging due to the doping dependent 
Berry phase associated with the JT effect |Q. How- 
ever, manganites at x = 0.5 exhibit a variety of ground 
states including the CE phase as in PrCa or NdSr, the 
A-phase as in PrSr or the metallic fcrromagnetism as in 
LaSr. A Monte-Carlo study of the two-band model with 
JT phonons ||2^ seem to capture most of these phases. 
An extension of our mean field theory incorporating the 



Jahn- Teller effect and breathing mode reproduces the 
charge/orbital ordered GE phase as well as the A- type 
phase [|0). 



IV. THE NATURE OF THE G-PHASE AND 
CANTING 

Expermentally ||^, Q it is seen that there seems to 
be little canting in the A and G phases. This was also 
emphasized by Maezono et. al. It is also seen 

that there is a predominant occupation of orbitals of 
one character in these phases. Recent experiments by 
Mahendiran et. al. on Smi-ajCa^MnOa suggest that 
even the G-phase for low doping may have little canting. 
The doped carriers seem to form ferromagnetic clusters 
leaving behind a uniform G-phase as background. In 
the band-picture, we have already noticed that for low 
electron doping the SE wins over the DE and the phase 
is G-type antiferromagnetic. One expects this phase to 
be canted as electrons gain kinetic energy due to the DE 
mechanism. The canting angle will be anisotropic, i.e., 
Oxy will be different from 0^ due to the anisotropy of 
the hopping integrals t^^p. However, no specific orbital 
ordering can be seen in this phase. This phase (without 
any orbital ordering) also has to be contrasted to the A- 
phase near x = 0.5 which has orbital ordering of dx2-y2 
type (see next section) . The stability of the G-phase near 
X = 1 is because of the dominance of antiferromagnetic 
energy whereas the stability of the A-phase near x = 0.5 
arises from the kinetic energy gain through DE in the 
plane due to selective d^2_y2 orbital ordering. Moreover, 
for finite Jh the canting is relatively small leading to a 
phase which closely resembles the G-phase at a; = 1. In 
Fig(5) we plot the canting angles as a function of Jn/t 
for a fixed value of JAp/t for some representative value 
of doping(x — 0.98). We find that the canting angle 
increases as a function of Jn/t for a given filling and 
JafI^ near x = 1. 

In the limit Jr 00, electron hopping to neighbour- 
ing sites with antiparallel core spins is not allowed. This 
is because the 'effective hopping parameter' for Jh — > 00 
is proportional to tcosiO /2) where 9 is the angle between 
the spins at neighbouring sites and antiparallel arrange- 
ment of spins reduces the 'effective hopping parameter' 
to zero. Hence the only way the electrons can take advan- 
tage of the kinetic energy gain due to increased doping 
is by canting the spins as much as possible. However, 
since ty's are anisotropic the canting angles will also 
be anisotropic. For a representative value of electron 
doping {x = 0.98) we find that there is no canting in 
the z direction and spins cant by about 10° in the x — y 
plane. This gives rise to a net ferromagnetic moment 
in the plane with a value higher than that across the 
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layers. Hence one would think of it as a canted A-phase 
as in iQ. However, inclusion of finite Jh changes this 
picture. A finite value of Jh allows the spins to go to 
^ wrong spin state^ at neighbouring site with an energy 
cost Jh- Hence the canting angle is reduced drastically 
compared to the Jh — *■ oo limit. In fact, for experimen- 
tally realistic values of Jh the canting is almost absent 
for low electron doping as can be inferred from Fig(5). 
(In fact one expects no canting for Jh ~ a,s DE is not 
operative.) Moreover, the kinetic energy gain which is 
proportional to the doping is also not effective in over- 
coming the SE energy. Hence one gets a canted G-phase 
with very small canting angles, thus resembling the G- 
phase at X = 1. Since the kinetic energy gain is also 
very small due to the smallness of the canting angle, this 
phase does not have any preferential orbital arrangement 
of the or d^2_y2 type as in the G and A phases. Thus 
we find that the stability of the G-phase is mainly due to 
the dominance of SE energy. This also means that the 
doping region over which the G-phase stabilizes will grow 
with increase in JAp/t- In particular, for JafI^ = the 
system should exhibit ferromagnetism for any doping 
making the G-C phase boundary collapse to the a; = 1 
point in the JAp/t-x phase plane. However JT4| find that 
the phase boundary between the canted G-phase and the 
C-phase does not change significantly as J afI^ is varied. 
More surprisingly, their phase diagram, if extrapolated 
to Jaf/^ = 0, will give the canted A-phase over a small 
region of doping near x = 1. In contrast to this, our 
phase diagram gives a ferromagnetic state for Jaf /t = 
for the whole doping regime and the stability region of 
the G-phase grows with increase in Jaf/^ in agreement 
with the physically expected result. Our results agree in 
general with the results of Maezono et. al though 
the A-phase near x ~ 0.5 is missing in that work. Sheng 
and Ting [p6| considered the problem from the strong 
correlation limit in contrast to the band limit which we 
have adopted. The C-phase between x — 0.6 and x = 0.9 
is missing in the strong correlation limit. 



V. ORBITAL STRUCTURE 

We find that in the C-phase the occupied orbitals are 
predominantly of character with a small admixture of 
dx2-y2. This happens because the electrons gain kinetic 
energy along the direction in which ferromagnetic corre- 
lations are stronger. For the same reason we find that in 
the A-phase the occupied orbitals are predominantly of 
dx2_y2 character. This, in effect, leads to suppression of 
hopping along antiferromagnetic bonds and explains why 
there is little canting in these systems. This is in agree- 
ment with experiments on Ndi_2;Sra;Mn03 Q. This also 
leads to a highly anisotropic band structure for G, C and 
A type structures and this feature becomes sharper as 



Jh increases. In particular, the G-phase has a quasi-one 
dimensional density of states. This also makes this phase 
very sensitive to substitutional disorder, possibly mak- 
ing it insulating. However, the A-phase is not sensitive 
to disorder and this rationalizes the (in-plane) metal- 
lic A-phase seen in experiments The nature of the 
occupied orbitals prevents electron motion along the z 
direction giving rise to a large anisotropy in the in-plane 
and out-of- plane resistivities. Experiments which probe 
the density of states, like tunneling measurements, will 
be able to see this feature. The low-temperature magnon 
spectrum will also throw light on the precise nature of 
the antiferromagnetic phase near a; = 1 and specifically 
the nature of canting in different manganites. 



VI. PHASE DIAGRAM OF THE ELECTRON 
DOPED BILAYER MANGANITES 

The present scheme of calculation can also be ap- 
plied to electron - doped bilayer manganites such as 
R2-2xAi+2xMn207 about which very little is known [^ . 
Since the interlayer coupling is roughly two orders of 
magnitude smaller than the coupling between bilayers 
one can apply the degenerate double - exchange, super - 
exchange model for a two layer system to study bilayered 
manganites. In this case the Brillouin zone is modified 
with kz taking only two values. As noted earlier the 
magnetic structure depends on the competition between 
the super-exchange and the kinetic energy renormalized 
by magnetic structure and orbital degrees of freedom. 
This suggests that in bilayer compounds where the ki- 
netic energy gain is predominantly in planes than in the 
z direction, the A-type antiferromagnetic phase is sta- 
bilized over the C-phase. This means that the dimen- 
sionality of the system plays a crucial role in the sta- 
bility of the C-phase. This can be clearly seen in the 
limit Jh oo where the band structure for G-phase be- 
comes one dimensional with e = — |icos(/c^). Detailed 
calculations support this picture as seen in Fig(6) where 
we present the results for a fixed value of JnS/t — 5 
and in Fig(7) where we present the results for a fixed 
value of JafS'^/I = 0.053. Battle et. al. [|7) have re- 
ported an A-type phase for NdSr2Mn207 [x — 0.5) and 
Ndi.iSri.9Mn207 {x = 0.45). We believe that this phase 
should extend even beyond x = 0.5 in accordance with 
our picture. Our phase diagram is in accordance with 
that of Maezono and Nagaosa (28) 



VII. DISCUSSION 

It is interesting to study the phase transitions between 
these anisotropic structures under an applied magnetic 
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field in z direction. We find that the G-type phase be- 
comes a canted A-type phase before transforming to the 
ferromagnetic phase for large x (close to 1). This is in 
agreement with recent experiments Further study is 
needed in this direction covering the whole doping regime 
0.5 < a; < 1. 

A major drawback of the current approach as well as 
that of earlier works is the homogeneous magnetic phases 
they predict. It seems likely that a phase separated 
regime is energetically more favourable than the canted 
phase 1^ . Phase separation, static or dynamic, seems to 
be a notable feature of manganites in the low-hole doped 
regime, charge ordered regime as well as the intermediate 
regime where there is a ferromagnetic metal to param- 
agnetic insulator transition as the temperature is varied. 
Batista et. al. through exact diagonalization studies 
of a single band model on small one dimensional clusters 
find that non-uniform ground states are highly possible 
in DE-SE systems. In particular, they find that at low 
electron doping, doped carriers get trapped at impurity 
sites and form ferromagnetic clusters. It will be interest- 
ing to study the two-band model to find exact nature of 
the G-phase near a; = 1. We expect the ferromagnetic 
clusters to be anisotropic in size with 'x — y radius' being 
larger than the 'z radius'. It should be possible to study 
phase separation using an orbitally degenerate version of 
the continuum model proposed by Soto et.al. It is 

also possible that the spiral |Q and the fiux phases ||3^] 
get stabilized for some values of doping as in the case 
of a single band double-exchange model though in our 
mean field picture we have not considered these phases. 
Work along these lines is in progress and will be reported 
elsewhere. 

To compare our results with the earlier work, we find a 
G-phase for low electron doping. We also find that the re- 
gion over which the G-phase is stabilized increases with 
JAp/t. This feature survives when canting is included 
as the canting angle is small for finite Ju/t- Our mean 
field theory takes into account the canting of core spins 
and also results in the A-phase near x = 0.5 (as seen 
in experiments) both of which arc missing in the work 
of Maezono et. al. jisj . Our model concentrates on the 
minimum number of relevant parameters and gives a uni- 
fied picture of the electron-doped manganites(including 
bilayers). This is in sharp contrast to the work of Mae- 
zono et. al. which uses five dimensionless parameters 
and separate order parameters for magnetic and orbital 
ordering. Our mean field theory also reproduces the C- 
phase between x = 0.6 and x — 0.9 which is missing in 
the strong coupling limit of Sheng and Ting ||2^ . We also 
clarified the nature of the G-phase near x = \ and the 
G-C phase boundary is as expected on physical grounds 
in contrast to van den Brink and Khomskii . 

In conclusion we have studied a model for electron- 
doped manganites with super-exchange between t2g elec- 



trons and double-exchange between orbitally degener- 
ate eg electrons. We find that finite Jh changes the 
phase diagram qualitatively. In particular the G-phase is 
favoured for low electron doping. This happens because 
the finite — Jh model, by allowing electrons to hop to 
neighbouring sites at an energy cost of Jh reduces the 
canting making the phase resemble more to the G-phase. 
The phase diagram agrees very well with the experimen- 
tal phase diagram of manganites for 0.5 < x < 1. By 
extending our mean field theory to incorporate the CE 
phase we find that it is stabilized over a wide range of 
values of JnSo/t and JApS^/t at x — 0.5. We extended 
this model for a two-layer system to predict the mag- 
netic phase diagram of electron doped bilayer mangan- 
ites. Here we find that the reduced dimensionality washes 
out the C-type phase. We also notice that the kinetic en- 
ergy gain due to DE leads to selective orbital ordering in 
the A and C phases while it is absent in the G-phase. We 
conclude that the present model qualitatively explains 
the anisotropic magnetic phases and believe that it can 
describe the phase transitions between these structures 
under an external field. A detailed study of this model is 
called for which should reveal the speculation about the 
phase separation in electron doped manganites. 
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Fig(l) : Phase diagram of the double-exchange and 
super-exchange model with degenerate Cg orbitals as- 
suming the doped electrons go into the F point and there 
is no canting of the core spins, (a) Phase diagram as a 
function of JafSq/I for a fixed value of JuSo/t = 5. (b) 
Phase diagram as a function of JuSo/t for a fixed value 
of Jaf5'^A = 0.053. 

Fig(2) : Phase diagram of the double-exchange and 
super-exchange model with degenerate Cg orbitals for a 
fixed value of JhSq/I — 5. Depending on the electron 
doping concentration and the ratio of the t2g superex- 
change to the Cg bandwidth, JafS'^/I, we find the A- 
type, C-type, G-type or ferromagnetic order. Values of 
Jhi t and Jaf were taken from density functional calcu- 
lation [9]. 

Fig(3) : Phase diagram of the double-exchange and 
super-exchange model with degenerate Cg orbitals for a 
fixed value of JafSq/I — 0.053. Depending on the elec- 
tron doping concentration and the ratio of the Hund's 
coupling to the Cg bandwidth, JHSa/t we find the A- 
type, C-type, G-type or ferromagnetic order. Values of 
Jh, t and Jaf were taken from the density functional 
calculation [9]. 

Fig(4) : Phase diagram of the double-exchange and 
super-exchange model with degenerate Cg orbitals at 
X = 0.5. 

Fig(5) : The angle difference between the neighbour- 
ing spins (in radians) for a representative value of doping 
X = 0.98 and Jaf/^ — 0.053 as a function of Jh /t- 

Fig(6) : Phase diagram of the bilayer system for a fixed 
value of JhSq/I = 5. Depending on the electron doping 
concentration and the ratio of the t2g superexchange to 
the Cg bandwidth, JafS'^/I, we find the A-type, G-type 
or ferromagnetic order. Note that the C-phase is missing 
in the bilayer system. 

Fig(7) : Phase diagram of the bilayer system for a fixed 
value of JAFSg/t = 0.053. Depending on the electron 
doping concentration and the ratio of the Hund's cou- 
pling to the Cg bandwidth, JhSq/I wc find the A-type or 
G-type phases. 
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